Entanglement within the Quantum Trajectory Description of Open Quantum Systems 
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The degree of entanglement in an open quantum system varies according to how information in 
the environment is read. A measure of this contextual entanglement is introduced based on quantum 
trajectory unravelings of the open system dynamics. It is used to characterize the entanglement in a 
driven quantum system of dimension 2 x cxa where the entanglement is induced by the environmental 
interaction. A detailed mechanism for the environment-induced entanglement is given. 
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Entanglement is a remarkable feature of quantum me- 
chanics that has attracted much attention in recent years 
for its potential use as a resource in quantum information 
processing 0. In practical situations a quantum system 
inevitably couples to its environment, so that the state 
of the system becomes mixed. Much effort has therefore 
been invested in finding a reasonable measure of entan- 
glement for mixed states. Specific proposals include the 
entanglement of formation |2, |^ and distillation , and 
the relative entropy of entanglement [3| . These measures 
are not practicably computable in general, however, and 
are readily accessible only for low-dimensional systems 
0, 01 and certain symmetric states . Recently, a vari- 
ational method for calculating the entanglement of for- 
mation (EOF) in general bipartite states was developed 
0. For bipartite mixed states, a computable measure of 
entanglement based on the negativity of the partial trans- 
pose has also been proposed 7] . In the case of continuous 
variable systems, the EOF of symmetric two-mode Gaus- 
sian states was recently obtained . 

In this Letter we consider mixed state entanglement 
from a new perspective. For pure states, the accepted 
measure of entanglement is settled. Considering a bi- 
partite state Pi = {ipi) {ipil (subsystems A and B), it is 
calculated as the von Neumann entropy 



E{pf) = -trB[pf\og2{pf)] , 



(1) 



where pf = trA(|V'i)(V'i|) '2]. It is well known that any 
mixed state p may be decomposed into an ensemble of 
pure states in an infinity of ways: p = TiiPilipi) {ipil, where 
Pi > and TiiPi = 1. Given a decomposition, the en- 
tanglement may be quantified by the ensemble average 
E = YjiPiE{pf). The perceived difficulty here is with the 
arbitrariness of the decomposition. Measures of mixed 
state entanglement differ in their strategy for replacing 
this arbitrariness by something specific. The EOF, for 
example, is defined as the minimum over all decompo- 
sitions; it intends to quantify the resources required to 
produce the state p under a specified quantum commu- 
nication protocol 2, 3]. The work reported in this Letter 
takes a different view. We observe that entangled states 
are pervasive in the quantum mechanics of composite sys- 
tems, often arising in situations that have no immediate 
quantum information connection. With the aim of under- 
standing entanglement in this broader sense, we propose 



to explore, rather than discard, the multiplicity of mixed 
state decompositions. The ambiguity arising in the mul- 
tiplicity is viewed as a signature of complementarity, i.e., 
an expected feature that refiects the very essence of a 
quantum mechanical description. 

We consider entanglement in bipartite open systems, 
specifically, where we resolve the arbitrariness in the de- 
composition of p by limiting attention to physically rele- 
vant decompositions — those for which pi is the probabil- 
ity of a classical record comprising information read from 
the environment. This information may be read in many 
ways. For each, there is a different decomposition of p 
and hence a different degree of entanglement. Quantum 
trajectory theory Q provides a natural measure of this 
entanglement in the form 



E,. 



(2) 



with p^.REC = trA(|V'c/;REc)(V'(7;REc|), where U labels a 
particular unraveling of the open system dynamics (read- 
ing of the environment) and RFC denotes a particular 
record; the overbar in Eq. (|2Jl denotes an average over 
records. The reading of the environment is continuous in 
time, and to preserve the purity of the state \ipU:REc), 
100% efficient (all scattered particles are ultimately de- 
tected). Ejj has a particular significance for each physi- 
cally realizable unraveling. 

We explore these ideas in an example, where we charac- 
terize the entanglement in a composite system consisting 
of (A) an optical cavity mode (harmonic oscillator), res- 
onantly coupled to {B) a two-state atom (single qubit). 
The atom is driven by a resonant external field and both 
subsystems are coupled to Markov reservoirs to account 
for the scattering of light into the vacuum of the electro- 
magnetic field. The Hamiltonian is 



H — Hab + 



(3) 



with 



Hab ~ fiuj{a) a + wb) + ihg{a)h — w a)^ 



iujt 



te*"*), 
- Har^ 



H 



BRe 



(4) 



where a and are cavity mode annihilation and creation 
operators ([a, d^] = 1), and b and are lowering and 
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raising operators for the atom ([&, = 1). The exter- 
nal field has amplitude 57, and A and B couple to reser- 
voirs Ra and Rb, respectively, via interactions Hara 
and Hbbs ■ The evolution of the reduced density opera- 
tor p in the interaction picture is governed by the master 
equation 



H^,„p] + {Ca+Cb)~p, 



(5) 



where Hab — iti{a''b — fa), H^^.^ = ifiil{b'^ — b), with 
ri = ^l/g, i = gt, and 



(6) 



where (o = a, b) and the damping constants Ta^b = ^a.b/g 
determine the reservoir interaction strengths. 

The dimension of A (g) i? is 2 x oo; hence the results 
of Refs. 0, IE IE 13 do not apply. The system exhibits 
the additional interesting feature of environment- assisted 
entanglement — i.e., an interaction with the environment 
is required to generate entanglement. Specifically, the 
degree of entanglement increases with the damping of 
the cavity mode A, passes through a peak, and then de- 
creases again. The behavior illustrates something of the 
subtlety associated with entanglement generation in open 
systems; the environment is not simply a source of de- 
coherence. Squeezing in the considered system shows a 
similar counterintuitive dependence on the cavity damp- 
ing Entanglement generation through an environ- 
mental interaction has also been reported for two cavity 
modes coupled to an incoherently driven atom [TH and 
for two atoms coupled to a common bath [T^ . 

We first uncover the origin of the environment-assisted 
entanglement. We note that Hab + H^^.^ has an eigen- 
state |q;)|(7), a — Cl, where \a) is a coherent state oi A and 
\g) denotes the ground state of B. The interpretation is 
that since the atom interacts with the sum of fields a — f2, 
it is stable in its ground state if the summed field ampli- 
tude is zero. In the presence of damping F;, 7^ 0, but with 
fa = 0, this eigenstate becomes the system steady state 
(a dark state) for all Ff, and f2 ^E|- Thus, if there is no in- 
teraction between subsystem A and its environment the 
system steady state is the product state and the 

entanglement Eu is zero; effectively the coupling between 
subsystems is turned off. The environmental interaction 
Fq destabilizes the dark state, restores the coupling, and 
generates entanglement. 

To describe the situation for F^ ^ it is convenient to 
introduce the displaced state p' , with p = (Cl) p' D {Cl) , 
where D{Cl) = exp[r2(a^ — a)\. Since _D(r2) is a local 
unitary operator, it does not change the degree of entan- 
glement. The displacement moves the driving from B to 
A 01 . In place of Eq. , we obtain 



dp 1 



H^^„p'] + iCA + CB)p' 



(7) 



where H^^^ = ihTa^ia - a'^). Note that F^ appears m 
two places: in the driving field interaction H^^, where 
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FIG. 1: Entanglement in direct detection, for = 1: typical 
temporal behavior of the von Neumann entropy for Fi, — 0.5 
and fa = 0.3 (a), 2.0 (b), and 20.0 (c); (d) entanglement 
[Eq. J^J] as a function of fa (for f i, = 0.5); (e) entanglement 
as a function of fa and f 



it governs the described creation of entanglement, and in 
the usual damping term Cap'- The interplay of the two, 
along with the interaction Hab , determines the behavior 
of the entanglement as a function of f q. In particular, in 
the limit F^ — > 00, the steady state of the cavity field in 
the displaced frame balances H^^ with Ca ■ This yields a 
coherent state of amplitude a' = — f2, the vacuum state in 
the original frame (see Fig. 01. Thus, in the large Ta limit, 
the steady state approaches p — paPb, where pA = |0)(0| 
and pb is determined by the balancing of H^^^ with C b , 
i.e., is the steady state of resonance fluorescence. In 
summary, for F^ = or 00 the entanglement is zero. 

Consider now the quantification of entanglement in be- 
tween these limits. Figure ^ illustrates the behavior with 
the unraveling of the density operator based upon direct 
detection of the scattered photons (unraveling U — D), 
those scattered through the cavity mirrors via the cou- 
pling of i? to A and those scattered directly by the atom 
B. Frames (a)-(c) illustrate the time dependence of the 
von Neumann entropy, where the entanglement is condi- 
tioned on a particular record of photon counts. Frames 
(d) and (e) show results after taking the average over 
records (a time average is permitted since trajectories 
are ergodic 0, 0). In the small and large damping 
limits, the entanglement is calculated to be 



-Alog2A-(I-A)log2(l-A), 



(8) 



with A « n*Tl and A « (2/Fa)2(4r2/Ffc)4, respectively. 
The peak in between these limits is shown in Figs, ^d) 
and (e), and is to be compared with the monotonic de- 
crease of Ejy with increasing f (Increasing f (, increases 
the rate of interruptions, due to the detection of a photon 
in the environment, that return the atom to the ground 
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state; each interruption destroys any entanglement cre- 
ated since the previous photon detection.) 

We turn now to the variation in the degree of entan- 
glement for different unravelings (U) of p. For simplicity, 
we set f fc = and focus attention on the behavior of 
the entanglement as a function of Ta ■ We begin by pre- 
senting an analytical approximation, which clarifies the 
behavior in Fig.^d) and serves as an introduction to the 
unraveling-dependence of the entanglement. 

Alsing and Carmichael 0| performed a semiclassical 
analysis of Eq. {7)). They found that as a function of the 
parameter ^ = {2Ta^)~^, the composite system exhibits 
a symmetry breaking transition at ^ = 1 (spontaneous 
dressed state polarization). If a = ^l{x + iy) and /3 = 
u-\- iv denote complex amplitudes of the cavity field and 
atomic polarization, in steady state 



(e>i), 
(e<i). 



(9) 



Figures|2Ia) and (b) illustrate the behavior of the steady- 
state amplitudes as a function of Ta- For Ta > l/2f2 
< 1), there are two permissible values for the ampli- 
tude of A, each correlated with an amplitude for B. With 
increasing Ta, the amplitudes of A separate and move in 
opposite directions on a circle of radius f2/2; they reach 
maximum separation at fa = l/\/2f2 and approach one 
another on the opposite side of the circle for Fa ^ 1. 
The correlated amplitudes of B move on a semicircle as 
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FIG. 2: Behavior of the semiclassical steady-state amplitudes 
a/Cl (a) and /3 (b) [Eq. @] with increasing Ta; the labeled 
points are (1) Ta = 0, (2) Ta = 1/2^, (3) Ta = l/^/2^, 
and (4) Fa S> 1. (c) Steady-state Wigner function for Fa — 
1 / \/2Cl (maximum peak separation) for = 3. 



shown. If f2 ^ 1, the phase-space separation across the 
circle is large compared with the vacuum-state uncer- 
tainty of oscillator A [point (3) in Fig. Efa)]. In such 
a case, the Wigner function for the reduced quantum- 
mechanical steady state pA = trB(p) is double peaked, 
as shown in Fig. Wi^) ■ 

The full quantum-mechanical steady state p is approx- 
imated, for f2 3> 1, by an equally weighted mixture of 
the states |q;_(-)|/3+) and |a_)|/3_), where \a±) are co- 
herent states of oscillator A and the states of B are 
\P±) = (\/^|e) + ^\g)) /V2, where \e{g)) denotes 
the excited (ground) state of the atom. Two compatible 
pure state ensembles are 



IV': 



REC, 



= e 



REC 



«+)|/3+) + (l-eREc)|a-)|/3-), (10) 



with 0REC = 0, 1 a dichotomous random variable (prob- 
abilities po — Pi ~ 1/2), and 



1^: 



REC) 



N{e 



"+>l^+> + e 



-i<S>T, 



(11) 



with 'I'rec uniformly distributed between and tt. For 
ensemble (|10|) the entanglement Ejj is zero; therefore the 
approximate p is separable. Nevertheless, the von Neu- 
mann entropy of state 1)111) is nonzero. We find it may 
be computed from Eq. ijHJl by replacing A with 



X ^ 1 1 VA^(1>R.Ec) + g^(<i>R,Ec) 

A±($REc) = 2±2 Bi^) ' ^ ^ 



where 



^(*REc) = Vl 



iin[2{n^xy + ^KEc)], (13) 



S($REc) =^ + e 



cos[2(f7\j/ + $REc)]. (14) 



After averaging over records a nonzero entanglement is 
obtained. It is plotted as the dashed curve in Fig.lSfb), 
where the behavior as a function of Fa resembles that 
in Fig. ^d), except Fig. |3Jb) shows a higher peak, ap- 
proaching unity for 1. Clearly there is a vast dif- 
ference in the entanglement computed for ensembles Hl()|l 
{Eu = 0) and ^ {Eu « 1). 

The difference has a physical origin; it arises from the 
way in which information is read from the environment. 
Ensemble (|1U|) yields Ejj — because it assumes that 
this information — encoded in 0rec — is able to distin- 
guish between \a+) and \a-). Ensemble lfTT|l assumes 
the opposite — that the record read in the environment 
cannot distinguish between |q!+) and |q;_). Direct detec- 
tion of scattered photons cannot distinguish these states, 
for example, because ~ \a^\'^ . Neither can mea- 

suring the amplitude Re(Q;-|-) — Re(Q!_). Measuring 
Im(a-(-) — — Im(a_), on the other hand, can. 

More generally, let us define the quadrature ampli- 
tude Xg — [e^^a^ + e~^^a) /2 and consider the homodyne 
detection unraveling of p (in the strong local oscillator 
limit; U — H) corresponding to a ^-quadrature measure- 
ment of the light scattered into reservoir Ra ■ The pure 
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FIG. 3: Entanglement in homodyne detection: (a) as a func- 
tion of Fa and Q, for 9 — (upper surface) and = Tr/2 (lower 
surface); (b) for fixed A = 3 and (i) = 0, (ii) 9 = 7r/40, 
(iii) 6 — it/10, and (iv) 9 — it/2, the dashed line is the entan- 
glement in ensemble IIH . All results are for Tt ~ 0. 



state ensemble is generated in this case by the stochastic 
Schrodinger equation 



yREC/ 



(15) 



where Hcs = Hab 
dqg = 



^cxt ^ Taa^a, and 



2ra{e"'a + e-'%^)TiEcdt + dW (16) 



is the record of charge deposited in the detector; dW is a 
real Wiener increment. Results for Eh are displayed in 
Fig. 131 In Fig. Efa) we plot the maximum and minimum 
values of Eh obtained with 9 = and tt/2, respectively, 
as a function of Cl and Fq. The minimum provides an 
upper bound on the EOF. It is not zero as in ensemble 
(fTn|l because Q is finite. In Fig. (|Sl)(b), the maximum Eh 
is well-approximated by Eqs. (|12(l - l|14|l even for fj = 3. 



The optimal entanglement occurs along a curve with fj 
inversely proportional to Ta, which is what one would 
expect from the definition of points (2) and (3) in Fig.|21 
Note, however, that with Si fixed, the peak, as a function 
of Fq, does not occur where |a+) and are maximally 
separated in phase space. This is explained by |/3+) and 
|/3_), which are not orthogonal at maximum separation, 
but only approach orthogonality as Fq ^ oo. 

Returning to our general theme, existing measures of 
entanglement view a mixed state p as a fundamental ob- 
ject and aim to associate a unique number with p. For 
open composite systems the mixing arises from a well- 
defined process. In this case, more can be said about 
entanglement by considering how the system acts upon 
its environment — nonlocally in the case of emitted light. 
Entanglement, as with the correlations it describes, be- 
comes a contextual notion. In this Letter we have shown 
how quantum trajectory theory can quantify such entan- 
glement, capturing the context-dependence in its differ- 
ent unravelings of p. 

Many questions are left open concerning observable 
ramifications of Ejj for different unravelings U. When the 
environment is monitored as a means of conditional state 
preparation lUj, physical implications are clear. More 
generally, do the records themselves contain indications 
that the system generating the measured outputs is de- 
scribed by an entangled state? The sensitivity of Ejj 
to minor changes in record making (imperfect detection) 
is also an interesting topic, with implications for deco- 
herence theory and the classical limit. These and other 
issues are left for future investigation. 
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